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Abstract 

We have investigated 3d electronic states of doped transition metals in II- 
VI diluted magnetic semiconductors, Zni-^^Mj^Se (M = Mn, Fe, Co) and 
Zni-j^Mn^iY (Y = Se, Te), using the transition-metal L2,3-edge X-ray absorp- 
tion spectroscopy (XAS) measurements. In order to explain the XAS spectra, 
we employed a tetragonal cluster model calculation, which includes not only 
the full ionic multiplet structure but also configuration interaction (CI). The 
results show that CI is essential to describe the experimental spectra ade- 
quately, indicating the strong hybridization between the transition metal 3d 
and the ligand p orbitals. In the study of Zni-^Mn^^Y (Y = Se, Te), we also 
found considerable spectral change in the Mn L2,3-edge XAS spectra for dif- 
ferent ligands, confirming the importance of the hybridization effects in these 

materials. 

PACS number(s): 71.23.An, 75.50.Pp, 78.70.Dm 
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I. INTRODUCTION 

There has been active studies on diluted magnetic semiconductors (DMS's) for the last 
two decades because of their unique magneto-optical and magneto-transport properties such 
as giant negative magnetoresistance, extremely large electronic g factor, large Faraday ro- 
tation, etclilfl Recently the interest has surged again because of their possible applications 
to "spintronics" based on the semiconductor technology.^ DMS's, which are made by sub- 
stitution of small amount of late 3d magnetic transition-metal (TM) atoms such as Mn, Fe, 
and Co ions into cation sites in the II- VI or III-V compound semiconductors, involve both 
semiconductor physics and magnetism. The band gap energy turns out to change with the 
doping concentration.Q Moreover, in some DMS's, the variation of the doping concentration 
was found to induce magnetic phase transitions, such as paramagnetic to spin glass phase 
and spin glass to antiferromagnetic phasejll and even ferromagnetism can also be induced 
by injecting carriers or photon irradiation.! 

Here our studies are mainly focused on the electronic structure of II- VI based DMS's. The 
II- VI compounds such as ZnSe and ZnTe have a zinc-blende crystal structure, in which the 
cation Zn^+ ions are under the tetrahedral T^ site symmetry. The TM doping induces merely 
small changes of the lattice-constants within the crystal structure. The (4sp)^ electrons 
of the TM atoms, which act like the Zn (4sp)^ electrons, participate in the ligand p-Zn 
Asp bonding-antibonding states. These states, which correspond to the respective valence 
band and conduction band of the host semiconductor, are hardly affected by the dopingElu. 
Meanwhile, the TM 3d electrons provide rather localized additional band states, which is 
expected to be located energetically in the wide band gap of the bonding-antibonding states. 
Optical absorption studies of DMS's showed the intra-atomic d-d transitions which can be 
understood by the 3d Coulomb multiplets excited from the corresponding ionic high-spin 
ground state.oQ In spite of the ionic characteristics, the TM 3d electrons are known to 
make strong covalent bonding with ligand p electrons. Furthermore, this covalent bonding, 
which is often represented by "p-(i" hybridization, is expected to play an important role 



in the interesting magneto-transport and magneto-optical properties such as giant Faraday 
rotation and Zeeman sphtting.El 

In order to understand the electronic structure of these DMS's, several photoemission 
spectroscopy measurements have been performed, mostly utilizing 3p -^ 3d resonance phe- 
nomena near the TM 3p absorption edge.tl't°l This technique enables us to separate out the 
TM 3d contribution to the electronic structure, and one can extract a sort of the 3d partial 
spectral weight. The resulting 3d spectral weight was found to be distributed in a very wide 
range of the valence band, no matter what the doping concentration is, indicating the strong 
p — d hybridization of ligand p band and the TM 3d orbitals. Meanwhile, the 3d spectral 
line shape does not agree with the band structure calculations.Ot3 For such reasons, the 
TM 3d partial spectral weights have been tried to be interpreted in terms of a many-body 
approach like a cluster model calculation with configuration interactionBily (CI) or an An- 
derson impurity model calculation.ll3 These model calculations, which are characterized by 
phenomenological physical quantities such as TM 3d on-site Coulomb energy, p — d hy- 
bridization, and the ligand p to TM 3d charge transfer energy, have been very successful in 
describing the electronic structure of TM compounds-M 

Soft x-ray absorption spectroscopy (XAS) is another powerful high-energy probe to in- 
vestigate the electronic structure of transition metal compounds.li^ Especially in diluted 
systems, XAS has big advantages over the photoemission measurements. Figure 1 shows 
XAS spectra at TM L2,3-edges of Zni_j.Mna.Se (x = 0.059), Zni.a^Fe^jSe {xd = 0.049), and 
Zni_a.Coa;Se {x = 0.050), which result from TM 2p — > 3d dipole transitions. The spectra 
are dominated by the large 2p core hole spin-orbit coupling energy, which divides them into 
roughly the L3 and L2 regions at low and high photon energy, respectively. The absorption 
process is a local process and its energy is determined by the characteristic core-hole energy. 
Despite the fact that the systems contain relatively small concentration of TM atoms, all 
the spectra nicely reveal the complex spectral structure, which is originated from the 2p 
core-hole-3rf Coulomb multiplets, and one can determine important physical quantities to 
understand the electronic structure by analysing them. 



The photoemission spectroscopy is known to directly probe the electronic structure and 
the TM 3p-edge resonant photoemission spectroscopy enables us to extract the TM 3d 
partial spectrum. However, the extracted spectrum somehow differs from the true TM 
3d electronic structure because of the different transition matrix elements in the resonant 
process from those in the direct photoemission process. Hence it is worthwhile to verify 
the physical quantities determined from the TM 3p-edge resonant photoemission studies 
by analysing the XAS results. Furthermore, in a view of the experimental technique, the 
XAS is rather bulk- sensitive, different from the resonant photoemission spectroscopy, and 
the results are much less sensitive to the surface condition of samples. Only a few XAS 
studies, however, have been performed for DMS's so far.ll3Lj Further, the reported studies 
were mostly focused on XAS spectra of the ligand atoms such as S K-edg&-n and Te Li-edge 
and L3-edge.lla The studies concluded qualitatively that the p-d hybridization is strong and 
the p-d hybridization strength depends on the ligand ions.ll3 However, quantitative analysis 
for the electronic structure has not been carried out. 

In this paper, we report high-resolution XAS studies to investigate the electronic struc- 
ture of various II- VI DMS's. The XAS spectra at the TM L2,3 edges are analysed by using 
a tetrahedral cluster model calculation including CI as well as the full ionic multiplets. 
Previously, the model calculations for TM L2,3-edge XAS spectra were performed only for 
condensed systems such as divalent Ni compounds,cy Ti203,Bll CoO,B3 and LiV02.E3, and 
showed that the CI considerably affects the XAS spectrum. Thus it is also important to 
elucidate how the effects of CI appear in the XAS spectra of diluted systems like DMS's. 

II. EXPERIMENT 

The samples used in this study were Zni-^^MuajSe {x = 0.059), Zni-j^Fe^jSe {x = 0.049), 
and Zni_2;Coa;Se {x = 0.050) thin films grown by molecular beam epitaxy on the GaAs (001) 
substrate and a Zni-ajMn^^Te {x = 0.60) single crystal with (110) surface. The details of 
the sample growth and characterization were reported elsewhere.a All the samples preserve 



the zinc-blende crystal structure. The XAS measurements were carried out at the Dragon 
high-resolution soft-x-ray beam lineo at the National Synchrotron Light Source, Brookhaven 
National Laboratory. The energy resolution of incoming photon was set to be about 0.5 eV 
in the full width at half maximum (FWHM). Before the measurements, the single-crystalline 
sample was cleaved in situ and the thin- film samples, which had been pre-etched using 1:3 
mixture of NH4OH (29 %):methanol,E3 were annealed in situ at 300 °C by radiation heating 
for about three hours.lH^ The measurements were performed at room temperature and the 
base pressure was maintained below 2x10^^° Torr. 

III. CALCULATIONAL METHOD 

A tetrahedral cluster model calculation, which we have employed to analyse the XAS 
spectra, includes not only the full multiplets of TM 3d electrons but also configuration 
interaction. The model is characterized by the on-site 3d Coulomb energy, U = E{d"'~^^) + 
E{d^^^) — 2E{d'^), the ligand-to-3rf charge-transfer energy, A = E{d^^^^L) — E{d'^), and the 
cation-d-ligand-p hybridization strength V. The initial configuration states are spanned over 
the ionic ground state 3d"' and the charge-transferred states, 3d"'~^^L, ■ ■ ■, 3d^^L}^^", where 
L denotes a ligand-p hole. The multiplets of the 3d electron configurations in the charge 
transferred states are taken into account explicitly while the total symmetry is preserved to 
be the same as the ionic ground state symmetry. The 3d Coulomb exchange interactions 
are represented by Racah parameters, B and C. Although the Racah parameters can vary 
slightly in different configurations, we imported the values from Ref. ^ and fixed them in 
all configurations for the sake of simplicity. The on-site Coulomb U and the charge transfer 
energy A are treated as control parameters in the calculation and defined with the lowest- 
energy 3d multiplet states of the corresponding configurations. 

Similarly, the final states are also spanned over 2p3d'^~^^, 2p3d"'~^'^L, ■ ■ ■, 2p3d^^I^~" 
configuration states. The Coulomb energies between the 2p core-hole and a 3d electron are 
represented by the 2p-3d Slater integrals, F^^, G^^, and G^^, and the spin-orbit coupling of 



the 2]5 core-hole is taken into account with (C2p)- The involved parameter values are borrowed 



from Ref. ^ and fixed for all configurations. The Slater integrals were scaled down by a 
reduction factor k in order to account for the solid state screenings.!^ The average Coulomb 
interaction Q between a TM 2p hole and a Sd electron is set to be an empirical value of 
1.2bU as widely adopted in the core- level spectroscopy studies. 

The hybridization interactions are expressed in terms of the Slater-Koster parameters 



Vpda and Vj,d.: Vt = {t\H\Lt) = ^^V^ + IV^d. and 14 = {\H\L,) = ^V^^n, where 
t{e) and Lt (L^) are TM-Sd and ligand-p orbitals with t2g (e^) symmetry orbitals under 
Td tetrahedral point group symmetry, respectively. For simplicity, we apply the empirical 
relation Vpda = —'^Vpdn-^ The crystal-field interaction lODq is fixed to be 0.25 eV in all 
the calculations, and the reason will be explained later. Here TM 3d spin-orbit coupling is 
neglected. It is known that the 3d spin-orbit coupling can be quenched when the ground- 
state orbital symmetry is either A or EiHa, and indeed, the ground state symmetries of Mn^+, 
Fe^^, and Co^"*" ions are ^Ai, ^E, and '^^2 under the T^ symmetry, respectively. 

The initial ground state is obtained in the modified Lanczos method. Then the XAS 
spectra are calculated by the continued fraction expansion.Ell The calculated spectrum is 
broadened with a Lorentzian function, which accounts for the core-hole lifetime, and fi- 
nally with a Gaussian function for the experimental resolution (0.5 eV). The Lorentzian 
broadenings, r2 for L2-edge and T^ for Ls-edge spectra are presented in Table I. 

IV. RESULTS AND DISCUSSION 

A. CI effects on the cluster model calculation 

The calculated TM L2,3-edge XAS spectra for Mn^+, Fe^"*", and Co^^ ions, which are in 
a tetrahedral cluster, are presented in Figure 2 in comparison with the corresponding exper- 
imental XAS spectra of Zni_^Mna.Se (x = 0.059), Zni.j-Fej^Se {x = 0.049), and Zni_a;Coa;Se 
{x = 0.050), respectively. Now the photon energy is presented relative to that of the corre- 



spending main L3 niultiplet peak with the highest intensity. In order to investigate the role 
of configuration interaction, the calculation has been performed for the orders of the charge 
transferred states. 

In the zero-th order calculations, which are denoted by LP_, the charge transferred states 
are neglected, i.e. the ionic calculations. Meanwhile in the m-th order calculations, denoted 
by L"^ , the charge transferred states are included up to the configuration states with m ligand 
holes, i.e. the calculations with the configuration states of 'idP'., 3(i""'"^L^, ■ ■ ■, 3(i""'"™'L™' for 
the initial states and the configuration states of 2p3c^"+^ 2p3d''+^l}, ■ ■ -, 2p3rf"+'"+^L'" for 
the final states. The spectra are normalized by the intensity of the L3 main peak. For a 
given TM ion, the same parameter set is employed for the L^, ■ ■ -, the L^ calculations. 

When we take into account only the L} charge transferred states, we find that the 
calculated spectrum changes considerably from the L^ spectra, especially in the cases of 
Fe^^ and Co^^ ions.E3 But the L^ spectra are nearly the same as the L spectra, indicating 
that the CI effects on the spectrum converge very rapidly. It is because the energies of 
the 3(i"'"'"^L^ and 2p3d"'~^^L^ configuration states already become much higher than those of 
3rf""^^L^ and the 2p3d'^^^L} states, due to the strong on-site Coulomb interaction. Thus we 
will consider, from now on, only the result of calculations which include configurations up 
to d^'^'^L^ and pd^^'^L^. Here the parameter values are chosen to give the L^ spectra, which 
agree well with the corresponding experimental ones. 

As can be seen in Fig. 2, the spectrum is strongly disturbed by CI, especially in high 
energy region of the main L3 multiplet peak. The peak separation becomes reduced and the 
sharp peak structure is smeared out with CI. Although similar trend might be reproduced 
just by increasing the crystal field splitting lODg value without including CljHa we beheve 
CI is essential for proper interpretation of the XAS spectra for the following reason. When 
we increase the value of IQDq without including CI, the pre-edge structure is well separated 
from the main peak and grows rapidly. At the same time the L2-edge spectrum becomes 
very complicated except for a Co^+ ion,E3 neither of which is observed in the experimental 
spectra. Hence, IQDq value should be small for Mn^"*" and Fe^^ ions but large for a Co^^ 
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ion if one fits the experimental spectra without CI. It is, however, hard to believe that the 
lODq value, which reflects the hybridization strength in the calculation, changes abruptly 
along this series. ISj 

In the comparison of the calculated spectra with the experimental ones, the agreement 
is rather good in overall. But one can still recognize some minor disagreement, which is 
expected to arise from the uncertainties in Slater integrals and spin-orbit coupling param- 
eters and the solid state effects which have been left out in the cluster model. Each final 
state multiplet may have different lifetime, and thus for the better fit, one should apply an 
appropriate core-hole-lifetime broadening for each multiplet state. The solid state effects 
seem to be rather severe in the spectrum of a Co^^ ion, and thus we employed the reduction 
factor K smaller than other TM ions. 

The parameter values for the best fit are tabulated in Table 1. The ligand-to-3d charge 
transfer energy A of Mn, Fe, and Co DMS's are 5.0, 3.0, and 2.5 eV, respectively. Here 
A is defined with respect to the center of the multiplet. The decreasing behavior of A 
value along the series agrees with the decrease of the difference in the electro-negativity 
between TM and ligand atoms. These values are similar to them obtained from the d-d 
optical-absorption study.la The hybridization parameter Vpda are determined to be 1.0, 0.9, 
and 0.8 eV for Mn, Fe, and Co DMS's, respectively. The small decrease of Ypda with the 
increasing atomic number is also consistent with the expected contraction of the 3d electron 
wave function, but these values are somewhat smaller than those obtained from the analysis 
of photoemssion spectraBll3 This is probably due to the small reduction of the effective 
hybridization strength in the XAS final states induced by the strong Coulomb potential by 
the localized core orbitals as pointed by Gunnarsson and Jepsen.E3 Similar behavior is also 
observed in the analysis of the Ni compounds using the Anderson impurity model.l23 



B. Hybridization effect 

In order to study the details of the hybridization effect, we compared the XAS spectra 
of Zn1_3.Mn3.Se and Zni_3;Mn3.Te where different hgands will give variations in the charge- 
transfer energy A and the hybridization interaction Vpda- As can be seen in Fig. 1, the 
Mn Ls-edge XAS spectrum shows a simple and distinguishable multiplet structure, and 
thus is very appropriate to study the variation of the hybridization effect. Figure 3 shows 
the experimental XAS spectra of Zno.94iMno.o59Se and Zno.4Mno.6Te at the Mn Ls-edge 
together with the corresponding calculated ones.E3 The experimental spectra show four peak 
structure, a highest intensity peak (denoted by A), two high energy peaks at about 1 eV 
(denoted by B) and 3.5 eV (denoted by C) above, and a small pre-edge structure (denoted 
by D), for both compounds. This structure is mainly originated from the Mn 2p^3d^ XAS 
final state multiplets of the Mn^"^ ion. 

The spectra display the common multiplet structure, but one can still recognize consid- 
erable difference in the relative peak energies and intensities for different ligand-ions. For 
examples, the relative energy of the peak B is lower for Te-ligand ions, while the width of 
the peak C is narrower for the Se-ligand ions. The best fits for the multiplet of these two 
compounds were obtained with the parameter values of A = 5.0 eV and Vpda = 1-00 eV for 
Zno.94iMno.o59Se and A = 3.8 eV and Vpda = 0.90 eV for Zno.4Mno.6Te. Other parameters 
such as the on-site Coulomb interaction U and the crystal-field interaction 10-Dg are fixed 
to be 8.0 eV and 0.25 eV, respectively, since they are expected to be hardly affected by the 
variation of the ligand ions. In the picture, A and Vpda mainly determine the hybridization 
strength between TM-rf and ligand-p orbitals. 

As can be seen in the figure, the calculated spectra show good over all agreement with 
the experimental ones. Chemical trends of the parameter values are consistent with the 
electro-negativity of ligand ions and also with the nearest neighbor distances between Mn 
and ligand ions. According to the Harrison's scheme,L3 the bond-length (d) dependence of 
Vpda is d~^'^, but the obtained values yield d~^'^ dependence using a simple central-force 
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modelEB for the Mn-anion bond length. Such a disagreement, which was also observed by 
Larson et ai, may be attributed to the chemical difference between Te and Se ions.tj The 
agreement in the region of the second satellite is not as good as the first satellite, especially 
in Zno.4Mno.6Te. We suspect that this is probably due to the magnetic interaction between 
Mn^+ ions in this Mn-rich compound, which can affect the spectral shape considerably.E3 

It is worthwhile to understand how the XAS spectrum is affected by the parameters A 
and Vpda- Fig- 4 shows the XAS spectra as a function of A for a fixed Vpda = 1.0 eV in the 
left panel and the spectra as a function of Vpda for a fixed A = 4.0 eV in the right panel. 
First, the calculation results show that the position of the peak B is closer to the peak A 
and its intensity becomes smaller when the hybridization effect increases by decrease of A 
or increase of Vpda- Second, the structure C in the experimental spectrum seems to be a 
single peak structure. However, it turns out to consist of several multiplet states. Third, the 
pre-edge structure D is mainly induced by the crystal-field splitting lODq. This structure is 
absent at lODq = but becomes distinguishable for large lODq-valne.^ In our parameter 
set with small lODq, the structure D is barely recognized as in the experimental results. In 
practice, it is difficult to determine the exact value of lODq from the fitting of XAS spectra, 
and we simply fixed lODq = 0.25 eV just enough to show the structure D. 

V. CONCLUSION 

XAS spectra of various DMS's at the TM atom L2.3 edges are presented and described by 
a tetrahedral cluster model including full multiplet structure and configuration interaction. 
Due to the strong hybridization between the TM-3d and ligand-p orbitals, inclusion of CI 
with reasonable physical parameters gives much better results than the model without CI. 
Thus it is important to incorporate CI for understanding of the XAS spectra of DMS's. 
Hybridization effects on XAS spectra are also investigated varying ligand ions of Mn DMS's, 
and the change of the spectral shape is well explained by the cluster model with CI. 
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TABLES 
TABLE L Parameters used to calculate TM L2,3XAS spectra of Zni_2:Mn2;Se (x = 0.059), 

Zni-a^Fe^Se (x = 0.049), and Zni-j^Co^iSe (x = 0.050). See text for detailed explanations. All the 

values are given in eV. 

A V,da U Q B C F^^ Gla Gl^ C2p K T2 Lg 

Mn2+ 5.0 1.00 8.0 10.0 0.119 0.412 6.321 4.606 2.618 6.846 0.8 0.45 0.35 

Fe2+ 3.0 0.90 3.6 4.5 0.131 0.484 6.793 5.004 2.844 8.200 0.8 0.70 0.45 

Co2+ 2.5 0.80 4.8 6.0 0.138 0.541 7.259 5.397 3.069 9.748 0.7 1.00 0.50 
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FIGURES 
FIG. 1. Experimental XAS spectra of Zni_^Mna;Se {x = 0.059), Zni-^^Fe^^Se (x = 0.049), and 

Zni_a;Coa;Se (x = 0.050) at the TM L2,3 edges. 

FIG. 2. Experimental (dots) and Calculated (lines) L2,3 XAS spectra of Mn^+, Fe^+, and Co^+ 
ions in a tetrahedral cluster depending on the number of configurations employed. L" denotes a 
configuration, which contains a L" term in the initial and final states, up to which was included 
in the calculations. See text for detailed information on employed parameter values. 

FIG. 3. Experimental (dots) and calculated (lines) XAS spectra of Zni_2:Mn2:Se (x = 0.059) 
and Zni_2:Mn^Te (x = 0.60) at the Mn L3 edge. The vertical line is a guide for the peak position. 
For parameter values of the calculations, refer to text. 

FIG. 4. Calculated -L3 XAS spectra of a Mn^+ ion in a tetrahedral cluster. Left panel: As a 
function of A. Vpda is fixed to 1.0 eV. Right panel: As a function of Vpda- ^ is fixed to 4.0 eV. 
The vertical line is a guide for the peak position. 
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